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BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present innovation relates to procedure for the detection of stress state, wherein ambulatory 
heart period is measured and the derived signal is segmented into physiological states. The term 
"stress state" means herein also its opposite, "relaxation state". 

2. Description of the Prior Art 

Heart period is among the most commonly used parameters in physiological monitoring. The 
wide use of heart period is related, on the one hand, to the availability of electrocardiograph 
(ECG) acquisition device for noninvasive monitoring and, on the other hand, to central role of 
heart period in the autonomic nervous system function and sensitivity to several physiological 
states and conditions. Heart period (or, its reciprocal heart rate) forms a basis for different types 
of analyses and may be defined as the series of intervals between consecutive QRS-waveforms 
in the ECG-signal. Another method of deriving information on the time distance between con- 
secutive heart beats is the detection of heart beat intervals from heart pulse signal. 

The fact that heart period is a complex product of several physiological mechanisms poses a 
challenge to the use of heart period in applied contexts. This is especially the case within ambu- 
latory measurement, that is, measurement that is performed within natural, free-living condition 
or field tests, outside of controlled laboratory environment and protocols. However, the multide- 
termined nature of the heart period may also potentate a derivation of additional physiological 
measures from the heart period signal by means of decomposing a series of heart periods into 
separate components that have a physiological interpretation. 

It is well known that both branches of the autonomic nervous system (ANS), the sympathetic 
nervous system (SNS) and parasympathetic nervous system (PNS) influence heart rate. It is 
conamonly known that the activity of the SNS and PNS produce, respectively, an increase and 
decrease in the heart rate level. It is therefore no suiprising that much of the work on assessing 
physiological functions and states using information on heart beat signal often addresses changes 
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in heart beat as stemming from the influence of SNS and/or PNS. Unfortunately, it is usually 
very difficult to determine precisely the effects of SNS and PNS on the heart rate, since it is of- 
ten not apparent which branch of the ANS determines changes in heart rate and, in addition to 
these mechanisms, there are several other mechanisms altering the level of heart rate both di- 
rectly £uid indirectly, many of which are not well-known. 

The prior art has documented several research lines attempting to use heart rate variability 
(HRV) to quantify more selectively the activity of SNS and PNS. It has been documented that, 
especially the power of the so-called high frequency (HF) component of the HRV in the fre- 
quency region of 0.15-0.50 Hz provides information on the level of parasympathetic outflow to 
the heart. Unfortunately, although it has been claimed in some instances that the so-called low- 
frequency (LF) component of the HRV in the frequency region of 0,04-0.15 reflects SNS activ- 
ity, the effects of SNS on HRV are rather unclear and it is known that several other mechanisms 
also influence HRV and especially the LF component, including PNS, hormonal responses, 
metabolic adjustments, and blood pressure control. Thus, increases and decreases in the level of 
heart rate and HRV may be due to several sources and therefore, it may be only possible to inter- 
pret changes in the heart rate level and heart rate variability as being indicative of the activity 
level of SNS and PNS during controlled situation and preferably with the aid of other measures. 

The concept of stress refers generally in physiological domain to a state of heightened level of 
physiological activity without inomediate or apparent requirements for such arousal. In this 
document, we use the state of stress to indicate a body balance wherein tiie overall cardiovascu- 
lar function, as indicated in, e.g., heart rate and cardiovascular output, is substantially higher than 
the level that is required by immediate physical metabolic requirements. The physiological state 
of stress may be due to different sources, such as, for example, physical load (e.g., posture), 
physical condition (e.g., fewer), mental stress, low level of resources (e.g., a burnout condition), 
or emotional arousal. 

It is of note that the detection of stress relates also closely the metabolic processes such as oxy- 
gen consumption and caloric consumption, since cardiovascular output indices alone would 
falsely indicate that metabolic requirements have increased during a state of stress. 

Feedback and information on personal stress state and more generally, resources would be very 
helpful for many individuals to monitor and manage their stress levels, to avoid a state of bum- 
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out, and generally to maintain and enhance health condition. Rest and relaxation are 
important features of stress management, as they help to reduce stress and further buffer and ac- 
cumulate resources against the onset and adverse effects of stress. 

It has been well-documented in the scientific literature that a state of stress is associated with 
heightened SNS influence to the heart and lowered or diminished PNS influence to the heart 
(e.g., Porges 1992). It is also known in the prior art that, at rest during steady conditions, relaxa- 
tion is shown as lowered level of cardiovascular activity and in specific, a decrease in the level 
of heart rate, and an increase in the magnitude of the HF component of HRV is often found to 
associate with state of increased relaxation. Some prior work has been documented to take ad- 
vantage of the role of HR and HRV in stress and relaxation related phenomena (US Patents 
4,832,038; 4,862,361; 5,891,044; 5,941,837; 6,104,947; 6,212,427; 6,358,201). 

Despite this correlational relationship, there has been not very much progress in the detection of 
stress-related physiological states on the basis of heart period signal. There has been some prior 
work on using information on heart period and HRV to classify user states, in particular in com- 
bination with other physiological measures such as skin temperature. The prior work based on 
ECG acquisition has been focused mostly on the determination of clinical condition with using 
specific autonomic nervous system tests and is therefore very limited in their application to char- 
acterize behavioral and physiological states in normal life in connection with, for example, am- 
bulatory measurement (US Patents 6,358,201, 5,299,199; 5,419,338; 6,390,986; 6,416,473). For 
example, the work presented by Childre et al (*201) applies heart rate variabiUty parameters in 
biofeedback context. The described invention is well-applicable to a controlled situation (e.g., in 
laboratory, relaxation training) but is clearly not applicable to ambulatory monitoring, wherein it 
would be crucial to differentiate physical and emotional sources of reactions and responses. Ac- 
cordingly, if used in ambulatory settings one should use a method of manually selecting time pe- 
riods of, e.g., relaxation training, to the analysis, as the described work does not include any 
means of separating different types of physical contexts from each other. 

There is also some documented work on the use of ECG and heart period derived measures to 
detect certain physiological conditions, wherein typically one or more parameters are monitored 
and a threshold limit is set to signal a change in state (US Patents 5,267,568; 6,126,595; 
6,358,201). These solutions are necessary linodted in the content of classifying states and suffer 
from the fact that the signal value of the heart period and HRV parameters is not always the 
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same but rather, typically varies in combination with physiological states. In other words, they 
do not account for the state- varying (conditional) relationships between heart period, HRV pa- 
rameters, and physiological states. 

There has been some work on the modeling of state-varying relationships in physiological sig- 
nals. However, the prior work is typically not related to the deteraMnation of stress, may involve 
heart rate measurement but require the use of two or more physiological measures (US Patents 
5,810,014; 5,846,206; 5,902,250; 5,921,937). As an example of the above the work of Davis et 
al. (*014) presents a general approach for a model-based identification of states according to 
multichannel measurement of raw bio-signals such as electrocardiographic, electromyographic, 
and electroencephalographic data. The presented system is geared towards the detection of spe- 
cific abnormal states from physiological waveforms by using a specific model fitting approach. 
It is obvious to one skilled in the art that the efficiency of the presented state modelling proce- 
dure is highly dependent on the availability of repetitive multichannel data, and thus does not 
apply for the in-depth analysis and decomposition of one signal, such as heart beat, to differenti- 
ate specific body states with characteristic dynamics. Furthermore, it should be clear that the de- 
scribed work is clearly not applicable to the analysis of heart beat signal to differentiate states of 
emotional stress fi'om other sources of body stress. 

It is thus clear that, firom the point of differentiating different physiological user states and in 
comparison to the acquisition of only one signal, these approaches require more effort on the 
measurement of physiological signals and are therefore susceptible to involving more material 
costs and more restricted user protocols. More importantly, the referred work does not include 
any contribution to the identification of stress and relaxation, wherein the occurrence of physical 
activity has not been able to take account in the context of using heart beat signal as the only sin- 
gle input. 

As indicated above, the major problem in the operationalization, measurement and monitoring of 
stress using information on cardiovascular function, such as acquisition of ECG and heart beat, 
would be the detection and differentiation of the sources of decreased and increased cardiac 
function. This is especially evident with increased cardiac activity (e.g., as shown by increased 
level of heart rate and decreased amplitude of HRV), which may result, for example, fi'om in- 
creased state of stress, increased state of physical activity, or postural changes. 
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OBJECTS AND SUMMARY OF THE INVENTION 

The object of the present invention is provide an automatic analysis tool for the purposes of pro- 
viding information on the physiological state of the user on the basis of ECG or heart period 
measurement. Any method of deriving information on the time difference between subsequent 
heart beats (e.g., different sensors) may be used to form the necessary input, heart beat signal, for 
this innovation. More specifically, the object of the invention is to provide a procedure for the 
differentiation of different physiological states and especially to provide a state-detection based 
method of deriving a measure of stress on the b£isis of heart period data. The characteristic fea- 
tures of the invention are disclosed in the accompanying claims. 

The invention is based on several computational steps wherein the order of the computations has 
some constraints. The steps may be characterized as follows: (1) initial transformations of ECG 
and/or heart period signals; (2) segmentation of the heart period data into stationary segments; 
(3) detection of the segments associated with other-than-stress related increases in cardiac activ- 
ity, including physical exercise, physical activity, recovery from physical activity, and postural 
change; (4) detection of segments characterized by relaxed state; (5) detection of segments con- 
taining a potential stress state; and (6) combining information obtained in steps 3-5 to provide an 
overall index of stress. The procedure may also contain an initial set-up of parameters, such as 
minimum heart rate, wherein some properties of the state detection system are either inputted or 
determined automatically. 

According to the present system, the detection of stress state is based on the following physio- 
logical assumptions: There is sympathetic dominance in relation to the parasympathetic domi- 
nance, as indicated in heart beat parameters, and there is no evidence of physical activity, exer- 
cise, movement or posture influenced cardiac reactivity. 

The invention may be described as an expert system that consists of a sequence of computations 
and inferences and provides a novel approach to the detection of physiological states and espe- 
cially stress on the basis of heart beat signal. The types and methods of giving user feedback are 
dependent on the purpose of the heart beat measuring equipment and software in connection with 
which the classification and state detection procedure has been applied. 
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The invention may be applied to and in association with devices such as heart rate 
monitors and other wearable and mobile computing devices, other types of equipments for 
physical monitoring and especially involving ECG detection, and software products suited for 
the analysis of heart period signal. The present invention is useful in combination of any physio- 
logical monitoring in the area of monitoring, enhancing and optinoizing physiological resources 
to better manage healthy lifestyle and wellbeing, sports training and fitness, and working capac- 
ity. It may be especially useful in the monitoring and providing information on state of stress, 
where it may be used in the long-term monitoring of the resources, i.e., accumulation and at- 
tenuation of stress and resources to deal with stress. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1. An overall view of the sequences in the procedure. 
Figure 2. Initial transformations of the ECG and heart period signal. 

Figure 3, An illustration of the procedure for the detection of exercise, recovery from exercise, 
physical activity, movement and postural changes. 

Figure 4. Procedure for the detection of relaxation, stress state, and combination of information 
to describe overall body resources 

Figure 5. Example of the detected periods of stress and relaxation states during a measurement 
period 

Figure 6. Example of accumulated resources during workday 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
INTRODUCTION 

The innovation is described here with the aid of an example implementation. It should be noted 
that the described system is not bound to any specific model or specifications, but rather, differ- 
ent alterations, forms, and improvements are possible and are in line with the spirit of the 
innovation. Thus, the following merely contains a description of the preferred embodiments of 
the innovation. 

The system involves initially an estimation, input, or use of previously measured parameter val- 
ues to characterize individual parameter values that influence the detection of states within the 
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system. This estimation of the values is optional for the system. The values may be 
estimated by using foraiulas based on empirical data and user-inputted values, such as age, 
weight, height, and sex. The values may be also based on history values based on empirical 
physiological measurements, or they may be inputted by the user or expert. Examples of values 
that may be used on fliis context are maximum and minimum levels of oxygen consumption, 
heart rate level, or HRV level. 

Figure 2 shows an overview of the initial transformations of the ECG and heart period signal. 
The ECG-signal is transformed into consecutive R-R intervals by using sequential QRS- 
complexes as markers of the beginning and end of a R-R interval. Also other components of the 
ECG signal may be used to detect the R-R interval. The consecutive R-R intervals are scaimed 
through an artifact detection filter to perform an initial correction of falsely detected, missed, and 
premature heart beats. Other potential methods of deriving information on the intervals between 
consecutive heart beats would be based on the detection of intervals from pulse signal on the ba- 
sis of different detection algorithms. It should be therefore noted, that in principle, the present 
innovation applies to all domains wherein temporal information is obtained on the timings be- 
tween consecutive heart beats. 

Consecutive artifact-corrected R-R intCTvals are transformed into equidistant heart period time 
series by using a weighted linear interpolation of the R-R intervals. For example, a time domain 
sampling rate of 5 Hz miay be used. 

Initial features that may be used in the segmentation of the physiological state of the user into 
stationary segments are computed prior to the segmentation. The features may include several 
different components, including the following: spectral power (ms^) in the LF and HF frequency 
regions, information on respiratory period and ventilation as derived, for example, from the heart 
period according to Kettunen & Saalasti (PCT/FI03/00426), variation and inconsistencies in res- 
piratory parameters, and oxygen consumption as derived from, for example, heart period and 
heart period derived respiratory period or ventilation. If available, also external measures of 
physical activity may be used as a component for segmenting the physiological processes into 
stationary epochs. All these measures are computed on a continuous basis with preset or dynami- 
cally alternating window sizes, thus potentiating the capability for real-time use. 
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A univariate or multivariate set of features is combined to segment the physiological 
processes into stationary segments. In other words, the target of the present procedure is to dif- 
ferentiate segments wherein the properties of the physiological system are similar and consistent 
within-the-segment, using the information obtained from the derived physiological features. An 
example of an algorithm to perform this is the so-called generalized likelihood ratio test, which 
basically performs the minimization of the variance within the segments and maxunization of the 
variance between the (consecutive) segments. This process is controlled by threshold parameters 
that determine the sensitivity of the segmentation process. The use of multivariate set of parame- 
ters describing user state is recommended to gain more stability and reliability to the segmenta- 
tion of data into stationary epochs. 

The generalized likelihood ratio test is classically applied in model selection but it is also 
adapted to segmentation algorithms, e.g. the one presented by Fancourt and Principe. In this em- 
bodiment the function used to estibtnate the signal in the segment is median. The error is calcu- 
lated from the mean absolute error between the median and the signal within segment. For closer 
details of the algorithm see the artiqle by Fancourt et al. It is clear that the segmentation may be 
implemented in alternative ways, e.g., with a method described in *The biomedical engineering"; 
Bronzino, Joseph D,CRC Press, Inc.,1995. 

Generally segments are defined from heart beat signal with a chosen rule (e.g. Francourt) for 
segmentation, and at least one segment describing a physiological state with elevated cardiac ac- 
tivity due to physical workload and/or increased metabolic rate is identified (e.g. eq. 1) ) and ex- 
cluded, if exists, and segments other than the excluded segments are detected for a potential 
stress state, which is identified using a predetermined rule for the heart beat signal (e.g. eqs. 3 
and 4). In a simple manner the segmentation uses so called moving window method with pre- 
defined window length. 

In the following the procedure for detecting physiological state is presented. Selected physio- 
logical parameters are analyzed from each stationary segment to form a basis for inference on 
states that include increased cardiac activity. The detection of physical activity (i.e., metabolism 
related) induced states of increased cardiac reactivity is done and those segments are excluded 
before the detection of potential states and intensities of relaxation and stress. Figure 3 presents a 
flowchart on the procedure of detecting physical activity related states. 
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Increased level of cardiac activity (e.g., increased heart rate level or decreased heart 
period level) may be due to the process of recovering from physical exercise or from any physi- 
cal activity. The state of recovering from exercise may be detected by, for example, applying the 
method described in Saalasti, Kettunen, Pulkkinen; patent applications FIN 20025038, PCT, to 
determine the level of recovery requirements in the body and classify segments involving recov- 
ery demands higher than a certain predetermined threshold as recovery state. 

The system requires information on the intensity of physical activity and exercise to differentiate 
exercise related effects on cardiac system from, e.g., stress related effects. The transformation of 
heart beat level into proportional intensity of physical activity is done by applying for example 
information on the relationship of heart rate level to oxygen consumption. It is a very practical 
way to determine oxygen consumption from heart rate and external input, e.g. accelerator and 
then using the derived value of oxygen consumption to determine whether physical load exists 
(not a stress state). This external input may comprise one or more sensors detecting movement, 
such as accelerometer or simply it can be manual input. 

It is possible also to use information on initial parameters, such as maximum and minimum lev- 
els of heart rate and oxygen or energy consumption and heart rate variability to increase the ac- 
curacy of determining the intensity of physical activity. However, it is clear that the direct trans- 
formation of heart level to exercise intensity does not produce optunal results for the purposes of 
differentiating other-than-exercise related increase in cardiac activity fix)m exercise induced re- 
activity. 

The estimate of the proportional intensity of physical activity may be enhanced by including in- 
formation on the respiratory period and ventilation, as obtained, e.g., from heart period. Further- 
more, given that exercise has also certain pattern of changes in respiratory activity (e.g., exercise 
is associated with increased respiration, whereas other-than-exercise related cardiac reactivity is 
not often associated with similar respiratory reactivity) and HRV (e.g., a certain level of HRV is 
typically associated with certain level of heart rate during exercise), and temporal length (i.e., 
exercise and intensive physical activity is necessarily bound to have a certain temporal length). It 
would be also possible to use external sensors to detect the occurrence and intensity of move- 
ment and physical activity with, for example, accelerometer or skin temperature sensors, loca- 
tion, distance, ventilation, or oxygen consumption, to support the determination of exercise in- 
tensity related changes in heart beat signal. 
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A certain threshold value is preset to determine whether a certain segment involves exercise or 
intensive physical exercise. Given that no other reasons related to respiration, HRV, or temporal 
length of increased heart rate provides evidence that a potential exercise induced increase in car- 
diac activity is not associated with intensive physical activity, segments containing a mean inten- 
sity of physical activity higher than, for example, 50% of exercise intensity may be classified as 
intensive exercise. 

Movements, start-up of a physical activity, and postural changes have all metabolic requirements 
and they also increase heart rate level. It is ttius necessary to differentiate these effects from 
other, non-metabolic factors influencing cardiac activity. It is known that PNS often controls in- 
crease in cardiac activity at the onset of movement and decrease in cardiac activity at the end of 
movement. It is also known that standing up, which poses a so-called orthostatic reaction, is as- 
sociated with increase in the spectral power in the LF frequency region of HRV in relationship to 
HF power, which usually decreases when standing up. Based on these assumptions, a covariance 
parameter according to Equation 1 has been derived to index movement related changes in car- 
diac activity. "E" denotes an average value- 
Equation 1 

^ LFpow^ 



movement = a'Cov{iog{HFpow))+ E 



HFpow^ 



An increase in the level of the described parameter indicates the prevalence of movement related 
increases in cardiac activity and thus may be used to differentiate states of movement, onset and 
offset of physical activity, and postural changes. When this exceeds preset threshold value, 
physical activity is detected. 

Information on the derivative of heart rate or heart rate variability is used as an indicator of 
movement induced cardiac reactivity. Relatively fast increase and decrease of heart rate level 
indicates physical load (not a stress state). 

The detection of light physical activity is similar to that of intensive physical activity but the 
threshold values are different. However, the detection of light physical activity differs from that 
of detecting intensive physical activity in that it is expected that metabolic demands are substan- 
tially lower and thus related changes in cardiac activity also different. . In other words, the co- 
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variance term introduced in Equation 1 is used to differentiate light physical activity related 
increases in metabolism from cardiac reactivity that is due to different sources. It is also here of 
note that, in a similar manner to the intensive physical activity, additional sensors detecting 
movement occurrence and intensity could be also used as a source of information to differentiate 
heart beat segments that include light physical activity. 

The detection of movement and reactivity due to changes in posture is performed by using the 
covariance term presented in Equation 1 with preset threshold values and a period surrounding 
the covariance. The covariance term may be smoothed with, for example, a Haiming window. 

In general, it is known that exercise, physical activity, and movement all have their typical time 
frames that have minimal temporal requirements for recovery. For this purposes, it may be pos- 
sible to use a frequency or time domain measure to provide information on the temporal proper- 
ties of the changes in cardiac activity, thus potentiating the comparison of those with those typi- 
cal in the case of different types of physical activity. 

It is important to notice that, for example, as the changes in the heart rates during physical activ- 
ity are associated with several different types of changes, these changes may be combined in an 
automatic decision function that may be, for example, deterministic, heuristic, or based on multi- 
variate methods such as fuzzy systems and neural network. 

In one embodiment, the information on whether the heart rate segment is or is not detected as 
having increased physical load such as movement or postural changes is used as an information 
in a model that determines oxygen consumption or energy consumption level on the basis of 
heart rate information. In its simplest form, the formula predicting oxygen consumption or en- 
ergy consumption level would be conditional to the information on the segment, so that a differ- 
ent formula or a correction factor would be used depending on whether the segment would be 
defined as having movement or other form of physical load. In another terms, given that stress 
and movement both alter cardiac activity and yet have different metabolic responses, the infor- 
mation obtained with the present innovation may be appUed to the context of evaluating meta- 
bolic processes. 

It is emphasized that in this embodiment both the oxygen and energy consumption may be util- 
ized in the detection process. Furthermore, a movement sensor may give further information and 
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provide improved detection for the segmentation of physical load from heart beat 
signal, and thus may improve both oxygen and energy consumption estimation. It should be clear 
there may be also other applications. 

Figure 4 shows an overall view on the detection of relaxation and stress. Relaxation index for the 
segment is determined by the combination of heart period and HF power as illustrated in Equa- 
tion 2. 

The LP- and HF power estimates may be calculated with, e.g., short-time Fourier transformation 
(a.k.a. Gabor transformation), providing a generalization of a stationary Fourier mto nonstation- 
ary time series analysis. Also Cohen-class of time-frequency distributions may be utilized, e.g., 
smoothed pseudo Wigner-Ville transformation, or time-scale presentations, e.g., Wavelet trans- 
formation. Furthermore, to improve system accuracy and behavior the signal may be altered with 
several digital signal processing steps, e.g., digital filtering, artifact detection and correction, de- 
trending, logarithmic transformations or signal scaling and normalization. 



The relaxation power may be detected, e.g., by calculating a mean HF power respect to heart rate 
inside the detected segment. A hard limit threshold may be defined to declare the reasonable 
function range for relaxation by exploiting the knowledge of the scale of the variables used in the 
calculus, e.g., a tiireshold proportional to minimum and maximum heart rate and HF power. An 
example relaxation power is presented in Equation 2. 

If no relaxation state is detected, the data segment is a potential candidate of including stress. 
Equation 3 shows an example of how an index of stress may be computed inside the segment. In 
Equation 3, HR denotes heart rate level, CT denotes inconsistencies in the frequency distribution 
of HRV due to changes in respiratory period, or altematively, variability in the respiratory signal. 
HFpow and LFpow denote spectral powers in the HF and LF regions of the HRV, respectively. 



Equation 2 



RLXpow = E( 



HFpow 



) 



HR 



Equation 3 
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It should be clear to anyone familiar with the field that the formula presented for computing 
stress index is only illustrative and may be formed in several methods. Both indices shown in 
Equations 2 and 3 may take different forms. Both indices may be also used to derive continuous 
parameters, which may be useful in some contexts of deriving overall, long-time information on 
stress and relaxation. 

The principle is to search for periods of increased heart rate with other markers of potential stress 
and related in general to heightened level of SNS activity and diminished level of PNS activity. 
The rate of recovery of heart rate and/or heart rate variability from state of increased cardiac ac- 
tivity may be also used as an indicator low level stress and high level resources. 

The indices of relaxation and stress may be combined taking advantage of the proportion of dif- 
ferent states and the intensity of relaxation and stress states. An example of such combination is 
presented in Equation 4. 

Equation 4 

T T 
Total ^resources = Cj ^-^-KLXpow — c^ ^-^-STRpow 

where ci and C2 are constants, T is total time of the measurement or a reference time, e.g., 24 
hours, Tr is time classified as relaxation, Ts is time classified as stress, RLXpow the intensity of 
relaxation, which may be scaled in absolute values or proportional to maximal value, for exam- 
ple, and STRpow is the intensity of stress state (see the computation of stress index in Equation 
3), which also may be scaled in absolute values or in proportion to the maximal value. The inten- 
sities may be calculated via mean or median values of the corresponding time series indices. The 
combined index may be especially useful in the comparison of different days. The constant ci 
and C2 may be used to scale the stress and relaxation intensities to each other. 

Figure 5 demonstrates a simplified example of the detected periods of stress and relaxation states 
dining a measurement period. Note that information is gained on both the intensity (denoted as 
fixed constants a-f in the figure) and length (e.g., duration) of relaxation and stress periods. 
The time periods that are not detected as either stress or relaxation are either excluded from the 
analysis as having exercise or other types of increased physical load, or are not excluded from 
the stress analysis but are not detected as containing either stress or relaxation state. 
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An example of using the equation 4 with data presented in Fig 5. In this example the level of 
RLXpow and STRpow are both pre-scaled between 0 and 1, that is, proportional to their maxi- 
mal and minimal values. The maximal and minimal values may be either inputted or detected 
from the data. If non-scaled values would be used the constants ci and C2 could have substan- 
tially larger or smaller values to account for, for example, logarithmic differences in the scales of 
STR pow and RLX pow. 



Total _ resources = 1 ~ • 0-50 " 2 ~^ ' 0.30 « 0.047 

where RU^ow = average(e, f) relaxation intensities and STRpow = average(a, b, c, d) stress 
intensities 

In principle, a segment is identified as stress if heart rate is constantly at a high level and it is ap- 
parent that it is caused by lowered level of PNS, higher level of SNS, and there is no evidence of 
physically determined heart rate reactivity. In fact, given tiiat stress is associated with a tonic low 
level of PNS and that high level of PNS potentiates reactivity, the absence of PNS mediated re- 
sponses is also an indicator of stiress and may be taken more formally into account. If a person is 
highly stressed, it may not be possible to find states of relaxation even during long periods. 

Feedback on the results of the present procedure may be presented to the user either on real-time 
application or off-line, after the measurement It is possible to give feedback on states, on relaxa- 
tion and stress components, and overall resources. The feedback may be given in many forms, 
for example, in graphics or as a single parameter that is easy to interpret. 

There are several methods of setting up die stress quantification procedure to account for accu- 
mulative and dynamic changes in stiress state. The accumulation of stress level (and, relaxation) 
may be determined as follows, 

Accumulated_stresst = Accumulated_stresst.i + f(Accumulated_stresst.i, Stress_levelt) 

The same kind of methods may be also used in the case of available resources, which may botii 
increase (i.e., more relaxation than stiress) and decrease (i.e., more stiress tiian relaxation). For 
this purposes, it may be useful to determine tiie change in resources according to, in general, fol- 
lowing. 
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Change_in_resourcest = f(stress_levelt, relaxationjevelt), wherein stressjevelt and 
relaxationjevelt may be determined by STRpowt and RLXpowt, respectively. 

In its simplest form, changes in resources may be determined as 
Change_in_resourcest = relaxationjlevelt - stress Jievelt, 

in particular if the stress and relaxation levels are scaled between 0 and 1. In this example of 
scaling, the decrease in resources due to stress would be at its maximum comparable to the in- 
crease due to resources at its maximum. 

The accumulation of resources due to momentary decreases and increases in resources may be 
determined as follows, 

Accumulated_j:esourcest = Accumulated_resourcest.i + f(Resourcest-i, Change_in_resourceSt). 

Fig. 6 demonstrates the product of a calculation of accumulated resources, wherein relaxation 
increases resources and stress state decreases accumulated resources. In this particular example 
stress level and relaxation level are scaled as proportional to their maximum. 

One inherent benefit of the accumulated index is in its capacity to describe the accumulated dy- 
namics in stress and resources usage. Another, very significant benefit of accumulated index is in 
its use to scale the accumulation of stress and resources usage according to a pre-determined 
scaling. 

It is possible to scale the relaxation and stress components to satisfy a pre-defined criteria in 
terms of resources usage. For example, the parameters can be set according to a scaling that it 
would take a period of 150 minutes to recover, at full relaxation level, firom stress state that 
would have been acquired, with full stress level, in only 75 minutes. Another pre-set criteria may 
influence the accumulation of stress level and accelerated use of resources, wherein history in- 
formation on accumulated stress state may further magnify the effects of present stress state to 
further accelerate the use of available resources. 
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Yet another possible scaling may be used to relate the accumulation and use of stress and/or 
resources to a scale that describes a probability that a particular accumulated stress or resources 
state would have healthy or unhealthy physical consequences. 

Implementations of the invention can be a computer software in a personal computer, a heart rate 
monitor (wrist top computer), ECG -monitoring or pulse monitoring equipment such as a cardiac 
pace maker and an ergometer (a stationary bicycle) or other fitness exercise equipment. Gener- 
ally an implementation consists of a processing unit, a terminal, software and at least one input 
device. 



